Mon. Not. R. Astron. Soc. 000, [THH] (2009) Printed 10 February 2010 (MN lAT^ style file v2.2) 



The long rapid decay phase of the extended emission from 
the short GRB 080503. 



^ F. Genet^'^*, N.R. Butler^ and J.Granot^ 

«^ ^ Center for Astrophysics Research, University of Hertfordshire, UK. 
Racah institute of Physics, Hebrew University of Jerusalem, Israel. 
^ Astronomy Department, University of California, Berkeley, CA, USA 



10 February 2010 



ABSTRACT 

The gamma-ray burst (GRB) 080503 was classified as a short GRB with extended 
emission (Perley et al. 2009). The origin of such extended emission (found in about a 
quarter of Swift short GRBs) is still unclear and may provide some clues to the identity 
of the elusive progenitors of short GRBs. The extended emission from GRB 080503 
is followed by a rapid decay phase (RDP) that is detected over an unusually large 
dynamical range (one decade in time and ~ 3.5 decades in flux), making it ideal for 
studying the nature of the extended emission from short GRBs. We model the broad 
envelope of extended emission and the subsequent RDP using a physical model for 
the prompt GRB emission and its high latitude emission tail (Genet & Granot 2009), 
in which the prompt emission (and its tail) is the sum of its individual pulses (and 
their tails). For GRB 080503, a single pulse fit is found to be unacceptable, even when 
ignoring short timescale variability. The RDP displays very strong spectral evolution 
and shows some evidence for the presence of two spectral components with different 
temporal behaviour, likely arising from distinct physical regions. A two pulse fit (a 
first pulse accounting for the gamma-ray extended emission and decay phase, and the 
second pulse accounting mostly for the X-ray decay phase) provides a much better 
(though not perfect) fit to the data. The shallow gamma-ray and steep hard X-ray 
decays are hard to account for simultaneously, and require the second pulse to deviate 
from the simplest version of the model we use. Therefore, while high latitude emission 
is a viable explanation for the RDP in GRB 080503, it does not pass our tests with 
flying colors, and it is quite plausible that another mechanism is at work here. Finally, 
we note that the properties of the RDP following the extended emission of short GRBs 
(keeping in mind the very small number of well studied cases so far) appear to have 
different properties than that following the prompt emission of long GRBs. However, 
a larger sample of short GRBs with extended emission is required before any strong 
conclusion can be drawn. 
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1 INTRODUCTION 

GRB080503 was detecte d by the Swift Burs t Alert Tele- 
scope (BAT) on 2008 may 3 (|Perlev et al.ll2009l l. Its prompt 
gamma-ray emission presents a short (~ 0.32 s) intense ini- 
tial spike followed by an extended emission lasting several 
minutes, for a total duration of Tqq ^ 232 s. The first short 
spike duration in the 15 — 150 keV band is 0.32it0.07 s and its 
peak flux is (1.2 ±0.2) x 10"'^ erg cm~^ s~^ The count rate 
hardness ratio between the 50 — 100 keV and the 25 — 50 keV 
of the initial spike is 1.2 it 0.3, consistent with other short 
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Swift bursts, but also with some long bursts. The fluence of 
the extended emission (measured between 5 s and 140 s) in 
the 15-150 keV band is (1.86±0.14) x 10"^ erg cm"^, about 
thirty times that of the initial spike, higher than for any 
other short Swift bursts, but with in the range of suc h ratio 
measured for BATSE short bursts. IPerlev et al.] (|2009l ) found 
a spectral lag between the 50 — 100 keV and the 25 — 50 keV 
bands consistent with zero. All this led them to associate 
GRB080503 with the "short" (|Kouveliotou et al.lll993l ) class. 

The Swift X-Ray Telescope (XRT) started observing 
GRB080503 about 82 s after the burst, detecting a bright 
early afterglow that rapidly decayed (a = 2 — 4 with oc 
t~^) to below the detection threshold during the first orbit, 
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which makes a record overall decline of ^ 6.5 decades, with 
a steep decay clearly observed for ^ 3.5 decades (see figure 6 
of Perley et al. 2009) . Chandra detections at ^ 3 days after 
the GRB indicate the presence of a separate X-ray after- 
glow component ^ 10^ times fainter than the peak emission 
from the X-ray tail of the prompt extended emission. The 
gamma-ray lightcurve of the extended emission (which can 
be found with the X-ray light-curve of the extended emis- 
sion in figure [1]) presents two relatively well defined spikes at 
^ 28 s and ^ 36 s, followed by some less defined variability 
with several narrow spikes (see also figure 1 of Perley et al. 
2009). The X-ray RDP lightcurves present a first break at 
~ 180 s seen in both soft and hard bands; later on, the hard 
X-ray band clearly shows a steepening (to a temporal slope 
of about —4.2) that is only marginally seen in the soft X- 
rays. The simultaneous gamma-ray decay is much shallower, 
with a temporal slope only about —2. 

The rapid decay phase (RDP) of GRB080503 is a 
smooth temporal and spectral continuation of the extended 
emission in the gamma-rays (Perley et al., 2009). This sug- 
gests that it is the tail of the prompt emission (O'Brien 
et al., 2006; Butler & Kocevski, 2007). The most popular 
model to explain the RDP following the prompt emission 
of long GRBs is High Latitude Emission (HLE; Kumar and 
Panaitescu 2000). In this model, after the prompt emission 
stops, photons from increasingly larger angles relative to the 
Une of sight still reach the observer, due to the curvature of 
the (assumed to be quasi-spherical) emitting surface. The 
blueshift of such photons decreases with the angle from the 
line of sight, and therefore with their arrival time. This re- 
sults in a simple relation between the temporal and spectral 



indexes of the flux, a — 2 + /3, where (x. t 



that 



holds at late times (when t — to ^ At, where to is the onset 
time of the pulse and At its width) for each pulse of the 
prompt emission. 

The RDP following the extended emission of the short 
GRB 080503, however, shows strong spectral evolution that 
is different than that typically seen in the RDP following 
the prompt emission of long GRBs. Whereas in the latter 
case the RDP usually shows a monotonic softening of the 
spectrum with time (see Zhang, Liang & Zhang 2007 and 
references therein), the RDP in GRB 080503 shows a similar 
softening in the XRT energy range together with a hardening 
in the spectral slope between the XRT and BAT energy 
ranges. This behavior suggests the presence of two distinct 
spectral components. 

Such extended emission is observed in about a quarter 
of Swift short bursts (Norris and Gehrels 2009; Norris et al. 
2009). It usually lasts tens of seconds, and has been shown to 
be always softer than the initial short sp ike, however show- 
ing n egligible lags as the initial spike (|Norris and Bonnel 
l2006l ). The same authors also found that the extended emis- 
sion is softer than the initial spike, making it closer to long 
bursts, which are softer than short bursts. However, the ini- 
tial spike of s hort burst with and with out EE show simi- 
lar properties ([Norris and GehrelsllioogI ). Many authors re- 
cently turned to environmental constraints to explain differ- 
ences between short GRBs with and without EE, but with 
results at best uncertain, sometime contradictory (Rhoads 
2008; Troja et al. 2008; Pong, Berger & Fox 2009; Sakamoto 
& Gehrels 2009; Nysewander, Fruchter k Pe'er 2009). There 
is as of now no strong explanation for the mechanism be- 



hind extended emission. It is therefore very important to test 
whether HLE can indeed explain the RDP of GRB080503, 
as this would allow some comparison with results for the 
RDP following the prompt emission without extended emis- 
sion (Willingale et al., 2009). This will stress similarities and 
differences between the extended emission and the prompt 
emission itself, ultimately leading to a better understanding 
of the origin of extended emission. 

Genet k Granot (2009; hereafter GG09) have devel- 
oped a simple yet physical and self- consistent model for the 
prompt and high latitude emission. The very large dynam- 
ical range over which the RDP of GRB080503 is observed 
allows us to test its behaviour at late times, which is rarely 
possible in other Swift bursts. In order to test whether the 
extended emission is coming from a mechanism similar to 
the prompt emission itself, we apply the GG09 model to 
GRB080503 data by folding it through the response matri- 
ces generated from the four energy bands of data (0.3-1.3 
keV and 1.3-10.0 keV for the XRT, 15-50 keV and 50-150 
keV for the BAT) . Section § [2] summarizes the theoretical 
model used to fit the burst lightcurves. Section §[3] describes 
how the data are reduced and fitted and the results of our 
fits, and in section § [H we discuss the impUcations of these 
results. 



2 THE THEORETICAL MODEL 

The prompt and extended emissions and their tail (the 
RDP) are taken to be the sum of its pulses and their tails, as 
described in GG09. Each pulse represents a single emission 
episode, assumed to come from an ultra-relativistic (7 ^ 1) 
thin spherical expanding shell, turns on at radius Rq and 
turns abruptly off at radius Rf = Ro-\- AR where AR is the 
width of the emission region. The emission is assumed to 
be uniform over the emitting shell and isotropic in its rest 
frame. The observed flux is calculated by integrating over 
the surface of equal arrival time of photons to the observer 
(following Granot 2005 and Granot et al. 2008). The char- 
acteristic observed times of a pulse are the ejection time of 
the emitting shell (Tej), the radial time at Rq (To) and at Rf 
{Tf = To(l + AR/Ro)); Tf is also essentially the temporal 
width of the pulse. For AR/Ro < a few, the peak of the 
pulse is observed at Tpeak = ^ej + Tf. We consider that the 
peak luminosity of the shell evolves as L'^, oc R^.ln the case 
of emission by synchrotron mechanism from electrons accel- 
erated in internal shocks and cooling fast, a — 1, but in the 
following it will be kept free as fixing its value does not sim- 
plify calculations and it allows to check for deviations from 
these assumptions. In this framework we also have the peak 
frequency of the uFu spectrum evolving as Up oc R~^. The 
emission spectrum is assumed to be the phenomenological 
Band function (Band et al. 1993). The number of photons 
N per unit photon energy E, area A and observed normal- 
ized time T = (T — Tej)/To of a single emission episode is 
then (where we also define ff = {Tf - Tej)/To): 



dN 



dEdAdT 



(^,T^1) = T- 



{T^ff) 



a + 2 



where a is kept free to vary as fixing its value does not 
simplify calculations and 

is the Band function with a normalization constant Bnorm, 
with high and low energy photon indexes of —1 — I3h and 
— 1 — I3i respectively, where A/3 = (3h— (^i and z = {E/Eo)f. 
One can remark that the observed spectrum is a Band func- 
tion spectrum, as the emitted one, but with its peak energy 
sweeping through the observed bands: the peak of E^dN/dE 

is ^peak = {l-/3l)Eo/f (where Eq = ^peak(T = 1)/(1-A)), 

which decreases linearly with the normalized time T. At 
times T ^ Tf is observed the high latitude emission of the 
shell. In this part of the lightcurve of the pulse, the spec- 
tral slope -d\og{E^dN/dtdEdA)/d\ogE will evolve from 
1 — Pi — TEohs/Eo (where ^obs is the observed photon en- 
ergy) to 1 — Ph if ^obs < Eq. This creates a break in the 
lightcurve, which becomes very smooth if the observations 
are done in an energy band of finite width. 



3 DATA REDUCTION AND FITTING 

We download the raw, unfiltered Swift BAT and XRT 
data for GRB 080503 from the Swift ArchiveG. Our reduction 
of these da ta to science q ualit y light curves and spectra are 
detailed in IPerlev etHI (2009). Here, we subdivide the X- 
ray and 7-ray light curves into much finer time bins to study 
the spectral evolution. We begin with the BAT 15-350 keV 
and XRT 0.3-10.0 keV light curves, with each bin containing 
sufficient counts to reach a signal-to- noise ratio (S/N) of 3 
or greater. We then generate respon se matrices for eac h time 
bin, using the tools summarized in Per lev et al.l (|2009 ). Our 
reduction accounts for the spacecraft slew for BAT as well 
as photon pileup for the XRT, among other effects. 

Next, we further subdivide the BAT and XRT data each 
into soft and hard energy channels: 15-50 keV and 50-150 
keV for BAT, 0.3-1.3 keV and 1.3-10.0 keV for the XRT. 
We rebin the data (and group the response matrices accord- 
ingly) so that each hard and soft channel still has S/N ^ 3. 
The soft and hard channel light curves are shown in Figure (TJ 
To fit the temporal/spectral model of GG09 to the 2-energy- 
channel BAT and XRT data as a function of time, we fold 
the model through the response matrices at each time bin 
and minimize the total • Fitting is accomplished using 
Markov Chain Monte Carlo (MCMC) through the python 
PyMC packag4!l. 

At the beginning of the RDP the flux is expected to be 
dominated by the last pulse. However, at later times other 
pulses are expected to become dominant, and the contribu- 
tion of each pulse to the flux at late times (t— to ^ At) scales 
as ^ FpeakT'^^^ when /3 is constant among the pulses: higher 
(with a larger Fpeak) or wider (with a larger T/) pulses will 
dominate. The width of the pulse having a larger power 
(2 + /3 ^ 4 — 5 for /3 ^ 2 — 3) in the relative contribu- 
tion to the flux, it will be the most important parameter to 
find the dominant pulse. This is not true if the high energy 
spectral slope is varying among pulses, but we assume that 

^ ftp:/ /legacy. gsfc.nasa.gov/swift/data 
^ http://code.google.eom/p/pymc 
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Table 1. One component model: best- fit parameters 



one component model 



parameter 


value 


A 


-0.67 ±0.18 




1.13 ±0.08 


tf [s] 


55.35 ±2.99 


tpeak [s] 


66.56 ±3.30 


dR/R 


8.33 ±1.43 


log Eo(tpeak) 


0.59 ±0.09 


Nh [cm-2] 


(8.7 ±2.2) X 1020 


a 


0.67 ±0.17 


log cro (BAT and XRT) 


-1.27 ±0.08 



Table 2. Two- component model: best- fit parameters 



two component model 



parameter 


1^* pulse 


2^d pulse 




(fixed) 


log norm. 


-4.69 ± 1.19 


-0.89 ±0.47 


/3h 


0.63 ±0.08 


2.13 ±0.33 


tf [s] 


37.23 ±2.45 


125.67 ± 10.37 


tpeak [s] 


42.12 ±2.24 


176.35 ±9.11 


a 


1.311 ±0.26 


-0.38 ±0.19 


log Eo 


1.0 ±0.3 


0.3 ±0.1 


dR/R 


6.67 ±0.39 


Nh [cm-2] 


(2.0 ±0.4) X 102i 



log (70 -1.63 ± 0.10 



such variations are small enough so that the contribution at 
late times of narrow pulses are still small compared to larger 
ones. Given the shape of the BAT lightcurve of GRB080503, 
it thus seems most natural to explain it by one or two broad 
pulse (s) superimposed by narrower ones that would account 
for the smaller timescale variability observed (such as the 
spikes at 28 s and 36 s) , but whose contribution to the RDP 
would be negligible: the ratio of their width being of the 
order of 10/60 ^ 0.16, the ratio of their contribution to the 
tail of the prompt emission would be about 0.16"^ ^ 7 10~^ 
or less. The results presented here will describe only the 
broad(s) pulse(s) fitting and results, as in it lay the most 
important physics. 

3.1 Single Pulse Model 

We begin by fitting a single pulse BAT emission model 
to study the expected late-time emission in the XRT bands. 
Initial temporal parameters (rise time t/ and peak time 
tpeak) are chosen to crudely reproduce an envelope of emis- 
sion containing the BAT light curve. We then allow these 
and the other model parameters to vary (Table [1]). 

We find that a single pulse model tends to yield a 
mediocre fit to the data. The fit requires a systematic er- 
ror term (in addition to the measured error) of 30% in order 
to yield = z/. This error is not driven by the BAT time- 
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102 10^ 

Time Since GRB 080503 [s] 




Time Since GRB 080503 [s] 

Figure 1. One (top) and two (bottom) component model fits 
to the early GRB 080503 data from the Swift BAT and XRT. 
The scaling of the counts data to flux units (y-axis) is approxi- 
mate; however, the model fitting appropriately convolves the time- 
dependent input spectrum through the instrumental response ma- 
trices. In the bottom panel, dotted and dashed lines show the two 
components separately, while the solid curves show the sum of 
the two components The best- fit parameters are shown tables Q] 
and\^ 

variability alone (although this contributes). There is clear 
difficulty in fitting the smooth light curve break in the X- 
ray bands and also potentially in fitting the high flux in the 
BAT after t ~ 150s. However, we note that the error bars 
on the late-time BAT flux are large. The limits plotted in 
Figure [T] (top) for the BAT hard channel are at the 2a level. 



3.2 Double Pulse Model 

Because the single pulse model has difficulty recon- 
structing the slow BAT decline as well as the rapid X-ray 
break, we attempt fitting a model with two pulses. As dis- 
played in Figure [1] (bottom) , the data are considerably bet- 
ter fit with two pulses, the parameters for which are given 
in Table [2] In the fitting, we tie several of the parameters 
from pulse B to those from pulse A. We also fix the low en- 
ergy Band model index f5i — 0, which is typical for GRBs. A 
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Figure 2. Spectra of GRB080503 at different time interval dur- 
ing the RDP when both gamma-ray and X-ray are observed: 
80 - 100 s (black), 100 - 150 s (dark blue), 150 - 290 s (red) and 
290 — 345 s (cyan). The dotted curve show for comparison a Band 
function spectrum with (3i = 0, (3^ = 1.3, and Fiy(Epf^) oc E"^]^, 
which is roughly expected from HLE. 



large value for that index from the single pulse model may 
further suggest the inapropriateness of that model. 

The data are fit by a soft early pulse (peaking at t = 42s 
with ^0 = 1.0 keV) and a softer later pulse (peaking at t = 
175s with Eq — 0.3 keV). The BAT emission is dominated by 
the early pulse, whereas the late pulse provides the temporal 
break observed in the XRT bands around 185 s, which in this 
case is the beginning of the high latitude decay of the second 
pulse. 

We have attempted fits of 3 or more pulses to the data 
as well. These do not appear to significantly improve the fit, 
which even for the two pulse model requires a 20% system- 
atic uncertainty. Apart from the possible slow BAT decay 
at t > 150s, this extra uncertainty appears to be needed to 
account for light curve variability, which is not treated in 
the model. 



4 DISCUSSION 

The main problem in fitting the data with a single pulse 
is the very steep decay of the hard X-ray band combined 
with the relatively shallow decay of the BAT bands. In order 
to find a reason for such a behaviour, we extracted spectra 
at different time intervals (80-100 s, 100-150 s, 150-290 s 
and 290 — 345 s) along the RDP lightcurve (when data in 
the four BAT and XRT bands are available) to probe the 
spectral evolution (see figure [21 where from earliest to lat- 
est spectra the colors are black, dark blue, red and cyan). 
Although the data are marginally consistent with a Band 
function with constant slopes and decreasing ^peak, there is 
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evidence for a departure from a Band spectrum to a spec- 
trum with a concave shape beginning at t > 150 s. Such spec- 
tral evolution would strongly suggest that at times > 150 s 
a second component contributes to the flux; this would nat- 
urally explain the distinct X-ray and gamma-ray behaviours 
of the RDP, accounting for the rapid X-ray decline and slow 
gamma- ray decay. 

The fit with two pulses gives much better results - as 
expected since two physical components seem to contribute 
to the flux. In this two pulse modelling, the first part of the 
emission (t < 65 s) is accounted for by a first pulse, and 
the second part, including the whole RDP, by the second 
component. The X-ray break at ~ 185 s is the start of the 
high latitude decay of the second pulse, and the smooth 
break along the rest of the RDP is due to the sweeping up of 
the peak energy of the spectrum across the observed energy 
bands (as explained in §[2]). Due to the reappearance at late 
time (t > 300 s) of the first component above the second 
one, the steep decay in the hard X-ray band is still not very 
well accounted for. However, when looking at the 4-points 
spectra, one can see that the high and low energy slopes seem 
to change (on top of the appearance of a second component): 
from the earliest to the latest spectrum, the low energy slope 
seems to decrease with time (the low energy slope of the 
vFi, spectrum becomes shallower with time), and the high 
energy slope also seems to steepen. This may explain why 
the introduction of a second component alone (with fixed 
Band function parameters) cannot reproduce both the steep 
hard X-ray and the slow gamma-ray decline, even with a 
very steep high energy photon index [f5h ~ 2.13): evolution 
of the Band function slopes seem required. 

The slow gamma-ray decay could alternatively be due to 
a series of spikes that are reasonably narrow (with At < t), 
unresolved (due to the wide late time BAT time bins) and 
hard (dominating at soft gamma-rays but with a small X-ray 
flux), so that we see only their smooth envelope, while their 
tail would hardly contribute to the observed X-ray decay. We 
do not directly model this, however, since it would introduce 
far too many new free parameter, and thus not provide a 
very stringent test for such a model. 

We had kept a as a free parameter since it did not com- 
pHcate the modeling and could test for deviations from its 
simplest form. Indeed, none of the values obtained for the 
fits are perfectly consistent with 1. Considering only the two 
pulse fit (the one pulse fit not being good enough to have any 
physical significance) , the values of a for the first and second 
pulses (pulse number indicated by a subscript number) are 
ai = 1.31=b0.26 and a2 = -0.38=b0.19. Since the value a = 1 
was obtained under the assumptions of electrons cooling fast 
by synchrotron emission in a coasting outfiow, this means 
that at least one of the above assumptions is not true in the 
case of GRB080503. The most natural assumption to relax 
is that the strength of the shocks in the outfiow is constant. 
Assuming that it varies with radius, namely parametriz- 
ing the relative uptstream to downstream four-velocity as 
Uud oc (while the Lorentz factor of the emitting shocked 
region remains constant), then a = a{u) and d = d{u) where 
L', oc and oc R^ . The dependence of a and d on u 
are summarized in table [3J where we have considered two 
options for the magnetic field in the emitting region: (i) a 
field strongly amplified at the shock that holds a constant 
fraction of the internal energy (equipartition) , and (ii) a 



Table 3. Dependence of the parameters a and d (with L' , oc R^ 

and Up oc R^) on u (with the strength of the shocks in the outfiow 
being parametrized by the relative uptstream to downstream four- 
velocity as Urmud OC R'^ ) , and values of u and d for the fitted 
value a2 = —0.38 zb 0.19 of the second pulse of the two pulses 
model. 



Equipartiton field 




w < 1 


w > 1 


a{u) 


l-2u 


l-2u 


d{u) 


5ii- 1 


3u-l 


u{a2) 


0.69 ±0.095 


0.69 ±0.095 


d[u{a2)] 


2.45 ± 0.475 


1.07 ±0.285 


Advected field 




w < 1 


w > 1 


a(u) 


1-u 


l-2u 


d{u) 


4u-l 


Su-1 


u(a2) 


1.38 ±0.19 


0.69 ±0.095 


d[u{a2)] 


4.52 ±0.76 


1.07 ±0.285 



pre-existing magnetic field advected from the central source 
(scaling as R~^ upstream) and merely compressed at the 
shock (advected field). The values of the parameters u and 
d corresponding to the value a2 = —0.38 ± 0.19 of the sec- 
ond pulse are given as well (since the value of a for the first 
pulse is close to its fiducial value a = 1). One can note that 
the value of a inferred form the fit implies a value of d > 1, 
meaning that Up rises with radius at least linearly, which is 
very different from the basic model where Vp oc R~^ . We 
stress, however, that the parameters a and d affect mainly 
the rising parts of the pulses and hardly affect the decaying 
parts. Since we do not have strong constraints on the time 
variation of the peak frequency during the rising parts of the 
pulses, parameter values similar to the ones we infer may be 
able to reasonably fit the data. 



5 CONCLUSION 

We have explored the extended emission and very long 
RDP of GRB080503 in the light of the realistic physically 
motivated model for the prompt and latitude emission from 
iGenet and GranotI (|2009l l . Neglecting narrow pulses that do 
not affect the high latitude emission, we fitted the broad un- 
derlying envelope of the extended emission and its RDP, ex- 
pecting a simple H LE behavior, motivated by the results of 
IPerlev et al.l (|2009l l . The first attempt, using a single pulse, 
failed to account for both the slow gamma-ray and very 
rapid hard X-ray declines. A close look at the spectral evo- 
lution during the RDP (at times > 150 s) showed evidence 
of a strong spectral evolution: the emergence of a second 
spectral component dominating the high energies is spotted, 
as well as evolution of the spectral slopes of the (assumed) 
Band function spectrum. This led us to fit the data with two 
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pulses, the first one accounting for tiie gamma-ray envelope 
and decay phase, the second one accounting for the X-ray 
decay phase. The fit thus obtained is far better - even if not 
perfect, since we are still not accounting for the smaller vari- 
ability. Adding a third pulse does not improve the fit, which 
suggests two pulses are enough to account for the data (as 
suggested by the spectral evolution). However, the fit we ob- 
tained is still far from perfect (the hard X-ray decay is still 
too steep at the latest times), the model being unable to 
account for the evolution of spectral slopes, and the fitting 
results showing discrepancies with its simplest form. This 
points out that the simplest form of GG09 model for HLE 
is not detailed enough to account for the complex RDP of 
GRB 080503. The mechanism at work here may thus still be 
HLE, but since this model has some difficulties and requires 
many degrees of freedom in order to produce a reasonable 
fit, this may suggest that another mechanism might be at 
the origin of the RDP in GRB 080503. 

A previous attempt at fitting GG09 model to a short 
pulse with extended e mission was done on GRB050724 by 
IWillingale et al.l (2009). They also could not obtain a good 
fit of the extended emission in both BAT and XRT band, 
the model not being able to account for a slow BAT de- 
cay and a steep XRT decay during the RDP - similarly 
to what is observed with GRB080503. Such feature may 
thus be a characteristic of RDP from extended emission, 
but more bur sts should be studied b efore any conclusion 
can be drawn. IWillingale et all (|2009l ) argue that HLE is a 
valid explanation for the RDP of GRBs without extended 
emission. If there is indeed a difference in the RDP of bursts 
with and without extended emission, then the properties of 
extended emission would be different from the prompt emis- 
sion, which might suggest they have a different physical ori- 
gin. GRB080503 may thus be the first stone paving the way 
to a better understanding of extended emissions in Gamma- 
ray bursts by studying their RDP. However, one should keep 
in mind that, on top of having been studied more carefully 
than other bursts, the RDP of GRB080503 is uncommonly 
long, meaning that maybe the range on which rapid decay is 
usually observed is not large enough to show deviations from 
the high latitude decay. In this case this should be taken as 
a caveat for drawing conclusions from model fitting to the 
rapid decay phase of GRBs. 



Granot J., Cohen- Tanugi J. k do Couto e Silva E., 2008, 
ApJ 677, 92 

Kouveliotou C., Meegan C.A., Fishman G.J., Bhat N.P., 
Briggs M.S., Koshut T.M., Paciesas W.S. and Pendleton 
G.N., 1993, ApJ 413, LlOl 
Kumar P. & Panaitescu A., 2000, ApJL 541, L51 
Norris J.P., Gehrels N. & Scargle J.D., 2009, ApJ submit- 
ted 

Norris J.T. k Bonnell J.T., 2006, ApJ 643, 266 

Norris J.P. k Gehrels N., 2008, GRBs 2007: Proceedings 

of the Santa Fe Conference, AIP 1000, 280 
Nysewander M., Fruchter A.S. k Pe'er, A., 2009, ApJ, 701, 

824 

O'Brien P.T., et al., 2006, ApJ 647, 1213 

Perley D.A., et al., 2009, ApJ 696, 1871 

Rhoads J.E., 2008, ApJ submitted (arXiv/0807.2642) 

Sakamoto T. k Gehrels N. 2009, Gamma- Ray Bursts: Sixth 

Huntsville Symposium, AIP 1133, 112 
Troja E., King A.R., O'Brien P.T., Lyons N., k Cusumono 

G., 2008, MNRAS, 385, LIO 
Willingale R., Genet F., Granot J k O'Brien P., MNRAS 

submitted 

Zhang B.B, Liang E.W. k Zhang B., 2007, ApJ 666, 1002 



AKNOWLEDGEMENT 

NRB is supported through the GLAST Fellowship Pro- 
gram (NASA Cooperative Agreement: NNG06DO90A). J.G. 
gratefully acknowledges a Royal Society Wolfson Research 
Merit Award. 



REFERENCES 

Band D., et al., 1993, ApJ 413, 281 

Butler N.R. k Kocevski D., 2007, ApJ 668, 400 

Fong W., Berger E. k Fox D.B., 2009, ApJ submitted 

(arXiv/0909.1804) 
Genet F. k Granot J., 2009, MNRAS 399, 1328 
Granot J., 2005, ApJ 631, 1022 



